
ABSTRACT 


The ways in which intense convective storms interact 
with their environment are considered for a number of specific 
severe storm situations. A physical model of subcloud wind 
fields and vertical wind profiles is developed to explain 
the often observed intensification of convective storms that 
move along or across thermal boundaries. A number of special, 
unusually dense, data sets are used to substantiate features 
of the model i. GOES imagery is used in conjunction with objec- 
tively analyzed surface wind data to develop a nowcast 
technique that might be used to identify specific storm cells 
likely to become tornadic . It is also shown that circulations 
associated with organized meso-a and meso-$ scale (250 to 
2500 km and 25 to 250 km respectively) storm complexes may, 
on occasion, strongly modify tropospheric thermodynamic 
patterns and flow fields. 
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INTERACTIONS BETWEEN CONVECTIVE STORMS 


AND THEIR ENVIRONMENT 

R. A. Maddox, L. R. Hoxit , and C. F. Chappell 


Mesometeorology Group, Atmospheric Physics 
and Chemistry Laboratory 
Environmental Research Laboratories, NOAA 
Boulder, Colorado 80303 


1. INTRODUCTION 

The evolution and life cycle of severe thunderstorms and 
convective mesosystems is determined by complex interactions 
among dynamic, kinematic and thermodynamic processes on scales 
ranging from those of small cloud particles to those of mid- 
latitude cyclones. Once deep convection develops in the 
atmosphere numerous reactions and interactions occur. Some 
responses act to suppress or dissipate the convection while 
others may act to intensify or perpetuate storms by imposing 
organization on the convective system or by triggering new 
developments. Most research has (for lack of comprehensive 
data sets) focused upon perturbations and interactions that 
occur within clouds or in the boundary layer beneath and near 
the clouds. 

The research effort reported here examined three specific 
types of thunderstorm/environment interactions. The first of 
these was the interaction of convective storms with shallow 



baroclinic zones such as warm fronts or thunderstorm produced, 
cold air outflow boundaries. Storms that move across or along 
such baroclinic zones are often observed to intensify and 
become tornadic (Miller, 1967). The meso-B scale kinematic 
and thermodynamic features that might explain this preferred 
intensification have been examined for a number of cases. 

Thunderstorms exist in and move through local mesoscale 
environments that may favor intensification because of specific 
spatial distributions of moisture, divergence, vorticity and 
conditional instability. It was hypothesized that these 
favorable meso-B scale regions might be detected and monitored 
utilizing conventional surface observations. Mappings of 
favorable fields or regions could then be superposed on 
concurrent satellite imagery to identify specific thunderstorm 
cells or complexes that were likely to intensify. These first 
two efforts considered down-scale interactions in which meso- 
scale environmental features interact with convective scale 
features in ways conducive to storm intensification. 

The third type of interactions considered were up-scale 
feedbacks in which organized convective complexes affect and 
modify their background environment. One obvious way in which 
this occurs is through continued, organized convective over- 
turning which modifies vertical temperature structure and 
consumes latent instability provided by the larger scale 
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environment (Fritsch, 1975a). The convective circulation is 
completed in regions of compensating mass subsidence. If this 
compensating subsidence occurs on organized, larger scales 
significant areas of warmer environmental temperatures result 
from the compressional heating. Mesoscale pressure and wind 
fields may be altered by these thermodynamic perturbations 
in ways which produce mesolow pressure systems (Hoxit et al . , 
1976), trigger new convection, act to organize and intensify 
existing convection, or even to modify large scale flow fields 
(Maddox, 1979a). A number of specific events have been studied 
in detail and the results are presented in subsequent sections. 

2. THUNDERSTORM/THERMAL BOUNDARY INTERACTIONS 

It has long been recognized that the area of intersection 
of a line of thunderstorms with a shallow baroclinic zone (e.g. 
warm front, another active squall-line, or a cool air outflow 
boundary produced by prior storms) was a favored location for 
the occurrence of severe, tornadic storms. Operational fore- 
casting papers by Magor (1959) and Miller (1967) emphasized 
the importance of this frequently observed characteristic of 
mesoscale convective systems. Kuhn et al . (1958) showed in 
composite, pre-tornado thermal charts that the favored location 
for tornado occurrence was northeast of the axis of a surface 
thermal ridge within a region of strong temperature gradient. 
Purdom's (1975 and 1976) recent analyses of satellite imagery 
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have provided "photographic" verification of the character- 
istics noted by earlier investigators and have also shown 
that similar increases in intensity occur even when individual 
thunderstorms move across or along such boundaries. 

Detailed studies of the thermodynamic and kinematic 
mechanisms responsible for these observed characteristics have 
not been undertaken. This has undoubtedly been because of a 
lack of mesoscale observations in the vicinity of such thermal 
boundaries. It has been assumed that ambient moisture con- 
vergence and relative vorticity must be maximized in these 
zones; and therefore, the most intense storms should be 
expected within these regions. 

Storm Data reports of tornadoes were used in conjunction 
with 3 -hourly surface charts to identify a number of thermal 
boundary tornado occurrences during 1977. (In addition, an 
important event that occurred in 1979 is also shown.) Examples 
of surface conditions preceding the tornadoes are shown for 
nine of these events in Figs, la through li . The GMT times 
of tornado occurrence from Storm Data are shown with the 
tornado tracks. 

Figure la illustrates a case early in the spring when a 
number of short-track tornadoes occurred in Oklahoma along a 
warm frontal boundary. These particular storms occurred along 
a very strong thermal gradient and were moving rapidly north- 
eastward. As the storms moved into deeper cold air they 
decreased in intensity. The situation that spawned the 
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destructive Wichita Falls, Texas, tornado is illustrated in 
Fig. lb. A number of intense tornadoes occurred along a thermal 
boundary that had been produced both by a large-scale warm front 
and by cool air outflow from storms to the north of the boundary. 

A case in May is shown in Fig. lc. Several strong tornadoes 
occurred in Missouri along the southern periphery of a thermal 
discontinuity that had been produced by thunderstorm activity. 
These storms occurred far to the east of a strong dryline 
and no tornado events were reported within the warm sector 
ahead of the dryline. Figure Id illustrates a case in June 
when a single , long-track intense tornado occurred in south- 
eastern Wisconsin. This storm tracked southeastward along 
the intense thermal gradient. The thermal boundary was likely 
produced by a combination of cool outflow from prior storms 
and by cooling from Lake Michigan. Once again no tornadoes 
were noted in the hot, moist air mass along and ahead of the 
surface cold front. 

Figure le depicts a July situation in which four intense, 
but short-track, tornadoes occurred in northwestern Wisconsin 
along a warm frontal boundary. These storms were moving almost 
normal to the isotherms and thus moved into an increasingly 
hostile convective environment. 
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Fig. la.’ Surface mesoanalysis for 2100 GMT, 2 March 1977. 

Surface features are shown along with 2 mb isobars 
and 4°F isotherms. The arrows are streamlines for 
the observed surface winds. Tornado tracks are 
indicated along with reported times of occurrence 
(GMT). 
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Fig. lb. Surface mesoanalysis for 2200 GMT, 10 April 1979. 

A thunderstorm cold air outflow boundary is depicted 
by the squall symbol with barbs added to indicate its 
direction of movement. The open barbs denote the 
dryline. Other details as in Fig. la. 


7 





•' ■ ' 1 -j 

Fig. lc. Surface mesoanalysis for 1800 GMT, 4 May 1977. 
Details as in Figs-, la and lb. 
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Fig. Id. Surface mesoanaly 
Details as in Fig 
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A late summer (August) case in southeastern Illinois is 
shown in Fig. If. The single, intense, long-track tornado was 
apparently spawned by a storm cell moving along a thermal 
boundary produced by earlier storm activity. This case is 
quite similar to that shown in Fig. Id. Figure lg illustrates 
another August case that occurred in central Minnesota. Two 
tornadoes moved along, and slightly across, a strong thermal 
gradient associated with a stationary east-west surface front. 
It appears that lake effect cooling from Lake Superior helped 
maintain the strength of the front . Tornadic storms were not 
reported in South Dakota and Nebraska ahead of the strong cold 
front. A fall situation (October) is shown in Fig. lh. Four 
tornadoes were reported in eastern Tennessee along an outflow 
boundary produced by earlier storm activity. Once again these 
tornadoes were far removed from the very strong surface cold 
front. Figure li depicts a winter (December) tornado situation 
along the Gulf Coast. In this case a single, strong tornado 
moved through the north suburbs of Houston. The storm was 
moving across a weak thermal gradient along a warm frontal 
boundary . 

These situations share many common features that are 
important whether they are considered relative to forecasting 
problems or to satellite data interpretation. The tornadoes 
tended to be intense (F2/F3 or greater, see Fujita, 1971) and 
destructive, but there were few individual occurrences. These 
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tornado events could thus be classified as small outbreaks - 
the type situation that may be harder to correctly forecast 
than a major (>10 tornadoes) outbreak (Galway, 1975). The 
storms that moved across the thermal boundary tended to produce 
intense, but short-track, tornadoes. Whereas, storms that 
moved along or parallel to such a thermal boundary produced 
intense, relatively long-track tornadoes. In all of the 
events shown, the thermal boundary tornadoes were the only 
tornado occurrences even though there were strong cold fronts 
and drylines pushing into the hot, moist air masses. This 
indicates that meso-a and meso-g scale conditions were probably 
not conducive to "family" type severe storm outbreaks. Tornadic 
storms were only able to develop in the meso-g scale thermal 
boundary zone. Therefore, the local environment along such 
boundaries must be considerably different than that of the 
hot, moist air mass. 

Observational and theoretical studies suggest that in the 
vicinity of thermal boundaries a local environment .conducive 
to severe storms is produced by horizontal variations in the 
boundary layer wind profile. Surface wind speeds over land 
are typically 50-60% of the geostrophic value and the flow is 
toward lower pressure with an isobaric crossing angle of about 
35 degrees. The wind veers with height and the speed increases 
until geostrophic or gradient balance is reached - usually 
1000 to 5000 feet above the surface. The surface crossing angle 
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Fig. le. Surface mesoanalysis for 0000 GMT, 31 July 1977. 
Details as in Fig. la. 
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1500 GMT 
21 AUG 1977 
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Fig. If. Surface mesoanalysis for 1500 GMT, 21 August 1977. 
Details as in Figs, la and lb. 
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Fig. lg. Surface mesoanalysis for 2100 GMT, 26 August 1977. 
Details as in Fig. la. 
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Fig. lh. Surface mesoanalysis for 2100 GMT, 1 October 1977. 
Details as in Figs, la and lb. 
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Fig. 


li. Surface mesoanalysis for 1200 GMT, 13 December 1977. 
Details as in Fig. la. 


16 




and the depth of the boundary layer are highly dependent upon 
atmospheric stability. 


A direct effect of strong low level thermal gradients is 
a rapidly varying pressure gradient with height. In the case 
of warm thermal advection the veering of the wind is increased 
while the reverse is true for cold advection. These modifi- 
cations of the basic Ekman wind profile have been documented 
in numerous studies (e.g. Blackadar, 1965; Lettau, 1967; and 
others). The change of the pressure gradient with height 
alters the vertical shear from that which would be present in 
barotropic conditions . Within a turbulent boundary layer the 
change in shear modifies vertical momentum transports further 
altering the boundary layer wind profile. 

Sheppard et al . (1952) were the first to suggest such 
modifications. Theoretical studies by Blackadar (1965), Cattle 
(1971) and Arya and Wyngaard (1975) all suggested that the 
angle between the observed surface wind and the surface geo- 
strophic wind is increased by low-level cold advection and 
decreased by warm advection. An extensive observational study 
by Hoxit (1974) produced similar results. Moreover, Hoxit ' s 
results indicated that the surface stress and ageostrophic wind 
component throughout the boundary layer were a function of the 
orientation and magnitude of the thermal wind vector relative 
to the surface geostrophic wind. The net result is that the 
average boundary layer ageostrophic wind component (directed 
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toward lower pressure) is greater in cases of cold thermal 
advection than in warm advection. 

A conceptual model of sub-cloud wind profiles near thermal 
boundaries has been developed using three case studies, time- 
series of upper-air soundings, instrumented tower data, and 
the results of the observational and theoretical boundary 
layer studies cited above. Figure 2 shows a schematic repre- 
sentation of wind profiles for the lowest kilometer in three 
different air masses. This is a typical meso-a pattern often 
associated with severe thunderstorms. At point A (within a 
warm, dry and well-mixed air mass) the wind veers slightly with 
height and the speed increases slowly. In the hot, moist and 
conditionally unstable air mass (point B) warm thermal advection 
and surface friction produce a wind profile that veers and 
increases in speed with height. At Point C ( in a cool, moist 
thunderstorm outflow region) the cold low-level thermal advection 
acts to decrease the veering of the winds. Winds at this point 
have maximum speeds near the surface and veer slowly with 
height until a transition occurs into the warmer air mass above. 
Within this cool, moist air mass cross-isobaric flow is often 
greater than it is within the hot, moist air mass. The Fig. 2 
cross-section of wind profiles indicates that the vertical 
wind profiles are modified in a manner which acts to maximize 
meso-g scale convergence and cyclonic vorticity (mean for the 
sub-cloud layer) within a narrow "mixing zone" between points 
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Fig. 2. Schematic- representation of boundary layer wind 

profiles within a typical severe thunderstorm producing 
surface pattern. Surface features and isotherms (°F) 
are indicated- Full wind barbs are 10 kt . 


B and C. Note that between points A and B the mean sub-cloud 
layer convergence and vorticity are actually smaller in mag- 
nitude than they are for the surface winds alone. 

The point C profile represents conditions that resulted 
from prior thunderstorm activity. The thunderstorm produced 
airmass will typically be meso-g in scale and therefore the 
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flow regime should not be balanced. However, the field should 
be evolving toward balanced flow and effects of cold air advection 
would be present in the vertical wind profiles. 

Since precipitation also may play a major role in setting 
up and maintaining thermal boundaries (especially warm fronts 
and thunderstorm outflows), evaporation in the vicinity of 
such boundaries may increase absolute low-level moisture values, 
and thus produce local areas of enhanced convective instability. 
This could be especially true if the thermal boundary were 
retreating slowly northeastward. In this situation the southerly 
flow in the hot, moist air mass would move across wet terrain 
as it approached the thermal "mixing zone". 

These hypothesized, systematic perturbations in planetary 
boundary layer wind profiles and moisture contents may, in 
part, explain why storms often reach maximum intensity and 
become tornadic as they approach thermal boundaries. Satellite 
imagery and satellite derived low-level wind fields might, on 
some occasions, be utilized to locate wind shear zones associated 
with thermal boundaries (see Maddox and Vonder Haar, 1979). 

a . The 27 May 1977 case 

NASA conducted an Atmospheric Variability Experiment 
(AVE) on 27 May 1977 that included a number of additional upper- 
air soundings over the central United States (see Hill and 
Turner, 1977 for details of the AVE’s). This day was chosen 
for study because serial soundings were taken at Oklahoma City 
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and instrumented tower data from NSSL were available. Figure 3 
shows a mesoanalysis of surface conditions at 1200 GMT, on 
27 May 1977, over Oklahoma and north Texas. Two weak outflow 
boundaries stretched across Oklahoma from northeast to south- 
west. These boundaries did not become active again during 
the afternoon and were not notable except for the special data 
taken within the cool air outflow region. 



Fig. 3. Surface mesoanalysis for 1200 GMT, 27 May 1977. 
Pressure analysis is of altimeter setting. 
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Time series of wind, temperature and wetbulb temperature 
data from the instrumented NSSL tower (WKY Television Systems 
481 m TV transmitting tower located 39 km north of NSSL and 
10 km northeast of downtown Oklahoma City) are presented in 
Fig. 4. The meteorological instrumentation and operation of 
the tower was explained by Carter (1970). Cool outflow air 
spread over the tower site between 1200 and 1300 GMT with an 
inversion persisting until about 1700 GMT. The vertical wind 
profiles from 1400 GMT on exhibit little veering with height 
and backing is evident within lower layers from 1330 to 1630 
GMT. Tremendous modifications of the thermodynamic structure 
occur from 1600 to 1700 GMT as surface heating modifies and 
destroys the cool air dome. Wetbulb temperatures rise rapidly 
and are quickly mixed through the tower layer (the region with 
wetbulb temperatures _> 16°C is shaded). 

A time series of the Oklahoma City AVE soundings is shown 
in Fig. 5. The shallow dome of cool outflow air extended up 
to about 900 mb (slightly deeper than the layer shown in the 
Fig. 4 time series). The wind profiles at 1400 and 1730 GMT 
are very similar to the hypothesized profile for point C on 
Fig. 2. A rapid change in the thermodynamic structure of the 
boundary layer between 1600 and 1700 GMT is also indicated by 
the Oklahoma City data. Of most interest is the rapid increase 
in moisture and 0 g values that rapidly mix upward to the 825 mb 
level. At Stephenville , Texas, (approximately 350 km south- 
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Fig. 4 


. Time series of observations (60 second averages at 
the hour and half hour) from the NSSL tower on 27 May 
1977. Solid lines are temperature ( C) ; dashed lines 
are wetbulb temperature (°C); and a full wind barb 
represents 5ms 


southwest of Oklahoma City) AVE soundings that were not affected 
by thunderstorm outflows were taken deep within the warm, moist 





airmass. These soundings indicated a more gradual diurnal 
modification of a deeper moist layer. At 2100 GMT the mixing 
ratio values at Oklahoma City exceeded the Stephenville values 
from the surface to 850 mb. Evaporation of surface moisture 
left by the morning thunderstorms may be partly responsible 
for the elevated mixing ratios in the Oklahoma City vicinity. 



Fig. 5. Time series of upper-air data taken at Oklahoma 

City on 27 May 1977. Dashed lines are temperature 
(°C); solid lines are 6 OO ; dash/X lines are mixing 
ratio; and a full wind Sarb represents 5 m s“l. 


b . The 6 May 1975 (Omaha tornado) case 


On 6 May 1975, a number of intense tornadoes occurred 
in eastern Nebraska along a very strong thermal gradient 
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associated with a warm front. This day was also a NASA AVE 
day; however , the experimental area was located south of 
Nebraska and no sequential soundings were available in the 
storm area. Hourly surface data, conventional soundings, radar 
data, and GOES imagery were utilized to study this event. 
Figures 6 and 7 depict mesoanalyses of surface conditions at 



Fig. 6. Surface mesoanalysis for 1800 GMT, 6 May 1975, 

Solid lines are altimeter setting ; dashed lines are 
4 U F isotherms; and regions with dewpoint temperature 
> 65°F are shaded. 
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1800 and 2100 GMT, on 6 May 1975. A Pacific front and occlu- 
sion was moving slowly eastward across eastern Kansas and 
Nebraska. A very slow moving warm front stretched from north- 
eastern Nebraska east-southeastward across northeastern 
Missouri. There was a strong temperature gradient north of 
the front that persistent cloudiness and thunderstorm activity 
in southern Iowa and northern Missouri helped maintain. Another 
area of active storms remained over central Missouri during the 
analysis period. 

The tornado events are shown on Fig. 7 with their reported 
times of occurrence. The storms were generally moving across 
the isotherms and all but one tornado had short tracks. A 
storm in northeastern Nebraska apparently moved northward 
almost parallel to the isotherms and produced a long-track 
tornado. Three of these tornadoes were F3 or greater in 
intensity. The storms fed upon a very narrow tongue of moist 
air that had been pulled up into eastern Nebraska. As the 
cells moved northward over cooler, drier air they weakened 
and were no longer tornadic. The general surface pattern is 
similar to that of the model in Fig. 2. The most favorable 
meso-B scale environment for severe storms should according to 
the model, lie along the warm front from northeastern Missouri, 
to extreme southeastern South Dakota, which is exactly where 
the most intense storms occurred. Indeed, radar film from 
Grand Island, Nebraska, and Des Moines, Iowa, indicated that 
the cell which spawned the Omaha tornado developed in the 
warm air sector south-southwest of Omaha. The storm became 
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tornadic as it moved across the "mixing zone" over the 
southwest suburbs of Omaha. 



Fig. 7. Surface mesoanalysis for 2100 GMT, 6 May 1975. 

Tornado events logged at the National Severe 
Storms Forecast Center are indicated (dark triangles) 
along with the GMT time of occurrence. Other 
details are as in Fig. 6. 
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Figures 8 and 9 present 850 mb and 500 mb analyses for 
0000 GMT, 7 May 1975 (approximately two hours after the Omaha 
tornado) . At 850 mb the Pacific front was very weak and is 
indicated to be dissipating. The surface warm front was shallow 
with the 850 mb position well to the north. Thus the storms 
might have remained severe as they moved north-northeastward 
were it not for the very narrow tongue of warm, moist air that 
was feeding them. Persistent storms in Missouri and Illinois 
had generated a mesohigh circulation that was detectable at 
850 mb. At 500 mb an intense cutoff low was moving slowly 
northward over Wyoming. The Pacific front and severe storm 
outbreak were associated with a weakening short-wave trough 
that was rotating around the cutoff. 

Time series of surface observations for Lincoln and Omaha, 
Nebraska, are shown in Figs. 10 and 11. (Lincoln is located 
approximately 90 km southwest of Omaha.) At Lincoln the warm 
front moved north of the station between 1700 and 1730 GMT. 
Thunderstorms developed in the warm air sector just west of 
Lincoln around 1800 GMT. These storms moved northward ahead 
of the Pacific front which passed Lincoln shortly after 2000 
GMT. At the time of the Omaha tornado (2140 GMT) Lincoln 
reported a temperature of 78°F with a dewpoint of only 45°F 
and west-southwest winds at 15 knots. Later in the evening 
the entire system took on the characteristics of a warm occluded 
system with the cool, moist surface air over Iowa pushing back 
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S westward. The cool air moved westward across Lincoln at 

I about 0200 GMT. 



Fig. 8. 850 mb analysis for 0000 GMT, 7 May 1975. 

Fronts, troughs, and outflow boundaries are 
indicated, as are 30 m height contours (solid 
lines) and 2°C isotherms. 
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The Omaha time series shows that the warm front remained 


south of Omaha throughout the period and that the Pacific front 
did not penetrate eastward to Omaha. The region just north of 
the warm front was characterized by strong and gusty east- 
southeasterly winds. These winds became very strong as the 
thunderstorm activity began. Between 2200 and 2300 GMT the 
surface winds reflect the passage of a small scale cyclonic 
circulation in which the tornado was apparently embedded. 
Although the surface pressure fluctuates violently between 2100 
and 2300 GMT, the surface temperature and dewpoint remain 
relatively constant. At the time of the tornadic storm Lincoln 
is located in about the relative position of point A and Omaha 
is at about point C as indicated on Fig. 2. Tornadic thunder- 
storms and severe hailstorms occurred along the strong thermal 
gradient northeast of the axis of hottest surface temperatures, 
c . The 24 April 1975 (Neosho tornado) case 

During the NASA AVE IV (see Fucik and Turner, 1975) 
experiment three to six hourly soundings were taken over much 
of the eastern 2/3rds of the United States during the period 
0000 GMT, 24 April through 1200 GMT, 25 April 1975. An intense 
tornado (F4) struck Neosho, Missouri, at about 0040 GMT on 
25 April 1975. (Neosho is located about 25 km south-southeast 
of Joplin, Missouri, and about 145 km east-northeast of Tulsa, 
Oklahoma.) Serial upper-air soundings were taken at Monett , 
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Missouri (Monett is 40 km east of Neosho), during AVE IV. 

The AVE data were used in conjunction with surface data, radar 
data, and GOES imagery to study this situation. 




Fig. 9. 500 mb analysis for 0000 GMT, '7 May 1975. 

Short-wave trough position is shown by the 
heavy dashed line with other details as in 
Fig. 8. 
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Fig. 10. Time series of surface observations for 
Lincoln, Nebraska, on 6 and 7 May 1975. 
Winds are in knots with a full barb = 10 kt 







Mesoscale surface analyses for 2200 GMT, 24 April and 
0000 GMT, 25 April 1975, are shown in Figs. 12 and 13. The 
surface situation was very complex with a number of mesoscale 
features present at 2200 GMT. The thunderstorm cool air 
outflow boundary that stretched across northern Arkansas into 
southeast Kansas and the mesohigh over southern Missouri had 
both been generated by early morning storm activity that 
moved from central Missouri into northern Arkansas before it 
dissipated. A strong temperature gradient was present along 
the southern and western peripheries of this outflow boundary. 

A mesolow pressure system was moving southeastward along the 
Oklahoma/Kansas border. Very warm surface temperatures were 
noted near the center of this mesolow and a tongue of very dry 
air had been pulled across central Oklahoma into its circulation. 
Very moist conditions were indicated by the band of 65°F plus 
dewpoints east and south of the mesolow. Dewpoints were high 
along the western edge of the outflow boundary. 

A mesohigh/low pressure couplet and another outflow 
boundary were located southeast of Kansas City. These features 
were associated with a severe thunderstorm that was producing 
large hail. A slow moving, large scale frontal system curved 
from the Texas panhandle northeastward across central Illinois. 
Thunderstorm cells were rapidly developing in Oklahoma in the 
warm, moist air mass southeast of Oklahoma City, and southeast 
and northeast of Tulsa (indicated on GOES photo but not shown). 
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Fig. 12. Surface mesoanalysis for 2200 GMT, 24 April 1975. 

Pressure analysis (solid lines) is of altimeter 
setting; isotherms (dashed lines) are shown for 
4°F intervals; and the region where dewpoints 
were > 65°F is shaded. 


By 0000 GMT the large scale front had pushed south over 
central and western Oklahoma pinching off the narrow dry 
tongue, and a large, severe hailstorm was in progress south 
and east of Oklahoma City. Two important mesolows were present, 
one in central Missouri and one between Joplin and Tulsa. A 
tornadic storm was in progress near Miami, Oklahoma, east- 
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northeast of one of the raesolow centers. This storm produced 
the Neosho tornado a few minutes later as it moved eastward 
across the thermal boundary zone. Neither tornado had a long 
track. The outflow boundary moving down from the north may 
also have been interacting with the tornadic cell by the time 
of the Neosho event. The two thermal boundary tornadoes were 
the only two reported during the afternoon and evening although 
there were many severe thunderstorms reported over eastern 
Oklahoma, southern Missouri and western Arkansas. Several 
funnels and tornadoes were reported during the late night and 
early morning hours in eastern Arkansas and Tennessee. 



Fig. 13. Surface mesoanalysis for 0000 GMT, 25 April 1975, 
with details as in Fig. 12. 
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Analyses for 850 and 500 mb are presented in Figs. 14 
and 15 for 2100 GMT, 24 April 1975, approximately three hours 
prior to the tornadoes. Radar echoes from the 2035 GMT radar 
summary chart have been superposed upon the analyses. At 850 mb 
a cyclonic circulation was located over south central Kansas 
and the dry-line and large scale front were detectable. The 
outflow air over southern Missouri and northern Arkansas was 
shallow with a south-southwesterly flow of warm, moist 850 mb 
air overriding the cool surface air. A stationary, thunderstorm 
generated thermal boundary stretched eastward from central 
Missouri. At 500 mb one weak short-wave trough was moving 
eastward over the Ohio Valley while an intensifying short-wave 
was moving across the Central Plains. A broad area of shower 
and thunderstorm activity was associated with the eastern trough 
and a number of storms were developing along and ahead of the 
western trough. 

Time series of the surface conditions at Ponca City, 
Oklahoma (located approximately 125 km west-northwest of Tulsa), 
Tulsa, Oklahoma, and Joplin, Missouri, are presented in Figs. 

16, 17, and 18. (At 2100 GMT Ponca City, Tulsa, and Joplin 
were located at about points A, B, and C on Fig. 2.) At Ponca 
City winds became westerly and the relative humidity decreased 
from 1500 to 2100 GMT while the station was within the tongue 
of dry air. The synoptic scale front pushed south of Ponca 
City about 2300 GMT and cooler, very moist air moved over the 
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station. At Tulsa very strong and gusty southwesterly winds 
were reported during the early morning hours. The thunderstorm 
outflow boundary appears to have moved just south of Tulsa 



Fig. 14. 850 mb analysis for 2100 GMT, 24 April 1975. 

Frontal boundaries and thunderstorm outflow 
boundaries are shown, as are 2°C isotherms (dashed), 
30 m height contours (solid), and 2035 GMT radar 
echoes (shaded). 
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Fig. 15. 500 mb analysis for 2100 GMT, 24 April 1975. Short- 

wave troughs are shown as heavy dashed lines with 
other details as in Fig. 14. 


and then southwesterly and increased in speed. The relative 
humidity decreased markedly between 2100 and 2300 GMT as the 
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station was under the influence of the dry tongue. The 
synoptic scale cold front did not pass Tulsa until after 0500 
GMT on the 25th. The Joplin time series shows the influence 
of several mesoscale features. During the early morning hours 
the jump in surface pressure and the abrupt drop in surface 
temperature occurred shortly after 1000 GMT as the thunderstorm 
outflow boundary moved past the station. The temperature 
remained quite cool until about 1600 GMT when the low cloud 
deck dissipated. After this time the surface temperature and 



Fig. 16. Surface observation time series for Ponca City, 
Oklahoma, on 24 and 25 April 1975. 
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Fig. 17. Surface observation time series for Tulsa, 
Oklahoma, on 24 and 25 April 1975. 


the dewpoint increased slowly as surface heating modified the 
cool air dome. Until 2300 GMT the surface wind slowly veered 
and decreased in speed as the mesohigh and attendant pressure 
gradient weakened. At 2300 GMT the temperature at Tulsa was 
more than 10°F warmer than the temperature at Joplin. After 
this time the winds shifted to the east as the outflow boundary 
from the north eased over the station and thunderstorms began 
to develop in the area. Between 0000 and 0200 GMT the pressure 
dropped rapidly and then increased slowly as the tornado cyclone 
passed to the south of the station. The temperature fell slowly 
during the period of pressure fluctuations and the wind speed 
never exceeded 10 kt. 
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Fig. 18. Surface observation time series for Joplin, 

Missouri, for 24 and 25 April 1975. 

The sequential soundings taken at Monett, Missouri, are 
shown in Fig. 19. The squall-line passed the station at about 
1130 GMT bringing cool and dry thunderstorm outflow air (low 0 ) 
which extended up to 875 mb. After 1500 GMT very dry (low 0 e > 
air advected northeastward over the cool surface dome. A layer 
between the cool surface air and the hot, dry air above main- 
tained mixing ratios of > lOg/kg through the period. The 
vertical wind profiles at 1547, 1804 and 2100 GMT are very 
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similar to the point C profile shown in Fig. 2. After 1900 
GMT the boundary layer air mass modifies very rapidly as 0 e 
increases and mixes quickly upward. The surface temperature 
increases slowly (1800 to 2315 change of 'v + 5°C) while 



Fig. 19, Upper-air time series for Monett, Missouri, on 24 April 
1975. Isotherms (dashed) are shown for 2°C intervals; 
9 g (K) is analyzed in solid lines for 4K intervals; 
mixing ratio lines (dash/X) are shown for selected 
values, with the region of mixing ratio >_ lOg/kg shaded 
and the region of < 3g/kg lightly hatched; and winds 
are in m s“l with a full barb equal to 5 m s"l. 
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at the surface jumps dramatically (1800 to 2315 change of 
'v + 22K). These changes are very similar to those shown for 
the 27 May Oklahoma case. The dramatic change in 0 g is, in 
large part, due to rapid increases in the local mixing ratio. 
Table 1 shows mean values of for the 950-900 mb layer for 
a number of stations. 




Table 

1 . Mean 

value of 0 g (K) 

for the 

950 - 900 mb 

layer 



Topeka- 

• Ft. Sill 

Stephenville 

Monett 

Little Rock 

Shreveport 



Kansas 

Oklahoma 

Texas 

Missouri 

Arkansas 

Louisiana. 

2100 

GMT 

325 

3 34* 

336 

332 

333 

336 

0000 

GMT 

322 

345 

339 

349 

337 

338 


*1900 GMT value 

The tremendous increase in 0 g at Monett may, in part, be due 
to the effects of very warm air moving across the damp terrain 
south-southwest of the station. 

3. ANALYSES OF SURFACE DATA AND CONCURRENT GOES IMAGERY 

A number of researchers have related surface thermodynamic 
and kinematic fields to the development of severe thunderstorms 
(e.g. Hudson, 1971; Newman, 1972; Sasaski, 1973; Doswell, 1977; 
and Barnes, 1978a). Other investigators have utilized satellite 
data sets to diagnose the severe thunderstorm environment. 

Peslen (1977) and Houghton arid Lee (197 7) computed kinematic 
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fields using winds derived from tracking clouds in satellite 
imagery. Negri et al . (1977) used satellite derived winds 
and satellite sounder precipitable water values to estimate 
moisture divergence. Adler and Fenn (1977) utilized IR cloud 
temperature to relate the rate of expansion of cold cloud tops 
to severe storm events. Purdom (1976) combined satellite data 
and surface observations to refine subjective mesoanalyses and 
Maddox et al . (1977) combined objective analyses of satellite 
derived winds and surface mixing ratio to estimate moisture 
divergence. It was surmised that important insights into storm 
severity might be gained by simply co-mapping gridded analyses 
of surface parameters and GOES imagery. 

Barnes (1978b) suggested that a tornado cyclone's ability 

to produce vortices in the friction layer depended upon ambient 

. -3 -1 

vorticity exceeding a threshold of 10 s on scales of about 
25 km. Although features on such small scales are not detect- 
able in routine surface data, it was hypothesized that analyses 
of meso-B scale ambient vorticity fields might indicate the 
likelihood of the vorticity threshold being reached at smaller 
scales. To investigate this hypothesis the vorticity equation 
was rewritten retaining only the terms (divergence and advection) 
that may be computed using surface winds: 

|| =--(£ + f)'V*"V - V *V(£ + f). (1) 
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The observed vorticity is (£ + f) and the time it would take 

-3 -1 

the mesoscale environment to spin up to 10 s , due to 

3 E 

divergence and advection (with assumed invariant), was 
called the Vorticity Time Parameter (VTP) , given by: 


VTP = 


10 3 - (g + f) 


-(£ + f) V • V - V • V(£ +f) 


s . 


( 2 ) 


Although the mesoscale environment (on scales resolved by 

-3 -1 

surface observations) certainly won't spin up to 10 s , . *. 

relatively smaller values of VTP (negative values are not 
allowed) may indicate that mesoscale ambient conditions are 
more favorable for the occurrence of intense tornadic storms. 

This hypothesis was investigated for the Omaha and Neosha 
tornado situations. The surface observations were extracted 
from WBAN-10 forms and analyzed using an objective analysis 
scheme developed by Barnes (1964 and 1973). The analysis 
response for the examples shown retained > 80% of the amplitude 
of waves longer than 200 km. (A detailed explanation of the 
analysis procedure is presented in Appendix I.) 

GOES images for 1800 and 2000 GMT, 6 May 1975, are shown 
in Figs. 20 and 21. Corresponding surface mesoanalyses and 
contoured values of the VTP have been mapped onto the satellite 
images. At 1800 GMT the favorable VTP region stretches along 
the frontal zone from southern Nebraska into central South 
Dakota. Several thunderstorms are developing in this zone 
from east central Nebraska into extreme southern South Dakota. 
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The VTP analysis indicates that these cells should be monitored 
for the next 1 to 2 hours. This is valuable "nowcast" infor- 
mation since strong storms are also indicated over Iowa, 
Missouri, Kansas, and Oklahoma. 



Fig. 20. GOES full resolution, visible image for 1800 GMT, 

6 May 1975. Surface features and VTP contours 
(values of 12X10^ s and less) have been superposed 
upon the image . 
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Fig. 21. GOES full resolution, visible image for 2000 GMT, 

6 May 1975. Details as in Fig. 20. 

At 2000 GMT the VTP pattern had split and showed two 
very favorable regions . One favorable area existed behind 
the surface front in a region where very dry air was being 
pulled into the low pressure circulation. No thunderstorms 
existed in this dry zone and the thunderstorms in southeast 
South Dakota were moving away from this favorable VTP region 
and into increasingly deeper low-level cold air. A second 
favorable VTP area was indicated along the Missouri River in 
eastern Nebraska. Several strong thunderstorms are indicated 
in the region and these would be the storms of concern for 
the next 1 to 2 hours. Indeed a check of the tornado locations 
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and times shown on Fig. 7 shows that the VTP analyses, when 
used with the GOES imagery , accurately specified the regions 
of tornadic activity. 

GOES images for 2100 and 2300 GMT, 24 April 1975, are 
presented in Figs. 22 and 23. At 2100 GMT there was a small 
region of favorable VTP associated with the mesolow circulation 
in southeast Kansas. There were no thunderstorms indicated in, 
or near to, this area. The location of the tongue of hot, dry 
air over central Oklahoma is quite apparent in the image. By 
2300 GMT the VTP values had decreased dramatically in the 
vicinity of the mesolow circulation. One large storm was 
almost centered within the favorable VTP region while two 
smaller cells were located on the peripheries of the region. 



Fig. 22. GOES full resolution, visible image for 2100 GMT, 
24 April 1975. Details as in Fig. 20. 
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The large cell spawned the Miami and Neosho tornadoes an hour 


and an hour and forty minutes later. Once again the combination 
of VTP and GOES imagery provides a valuable nowcast tool since 
the visual appearances of storms over southern Oklahoma and 
northern Missouri were more impressive than that of the tornadic 
storm. 



Fig. 23. 


GOES full resolution, visible image for 2300 GMT, 
24 April 1975. Details as in Fig. 20. 
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Although only two cases have been examined, the results 
are encouraging and indicate that fields derived from conven- 
tional observations might be used in conjunction with satellite 
data to specify regions where the threat of intense tornadic 
storms is maximum. The nowcast forecasting and warning appli- 
cations could be quite valuable once such a technique is tested 
and refined for a variety of tornado situations. An interactive 
data processing and analysis system, such as AOIPS, would be 
necessary to compute, analyze, contour, and display these types 
of data blends in a real time mode. 

4. CONVECTIVE MODIFICATIONS OF THE ENVIRONMENT 

The preceeding two sections have considered ways in which 
the near environment acts to control the intensity and evolution 
of convective storms. This section presents results documenting 
the life cycle of mesoscale pressure systems and relates these 
features to thermodynamic feedbacks that were produced by 
convective storm complexes. Many of the findings are from a 
detailed case study of a major severe thunderstorm outbreak 
that utilized satellite, radar, surface, and conventional 
upper-air data along with a set of special 1800 GMT sounding 
data . 

a . The 4 April 1977 (Birmingham tornado) case 

On 4 April 1977, a major outbreak of severe weather 
occurred over the southeastern United States. Twenty-two 
tornadoes were reported including one of F5 intensity at 


51 



Birmingham, Alabama, that killed 22 people while inflicting 
tremendous damage. An intense storm, with associated large 
hail and tornadic activity, was partly responsible for the 
crash of Southern Airways Flight 242 in northern Georgia and 
the deaths of 71 persons. Special 1800 GMT rawinsonde data 
were taken over the southeast U.S. as part of a SESAME ’77 
experiment (see Kreitzberg, 1977). These soundings were used 
in conjunction with surface observations, radar data, and GOES 
imagery to prepare detailed analyses of the meteorological 
conditions associated with the event. Miller (1978) has 
also analyzed this particular severe storm situation. 

Mesoanalyses of surface conditions over the Southeast are 
presented for the period 1200 GMT, 4 April to 0000 GMT, 5 April 
1977, along with selected satellite images, in Figs. 24a 
through 24i. Analyses of mesoscale pressure perturbations 
(large scale pressure subjectively estimated from time series 
data and then subtracted from the observed pressure) are on 
the 1200, 1500, 1800, 2100 and 0000 GMT charts. At 1200 GMT 
(Fig. 24a) considerable thunderstorm activity was already in 
progress, especially over Tennessee and Mississippi. A 
mesolow pressure center was located over northwestern Mississippi 
and a squall-line trailed southwestward into eastern Louisiana. 
The thunderstorm activity was occurring well to the east of the 
synoptic scale cold front that was just pushing into western 
Arkansas. Moist (dewpoints _> 65°F) southerly flow was indicated 


52 



from the Gulf Coast northward to Tennessee. As the day 
progressed the thunderstorm outflow boundary over western 
Mississippi moved slowly eastward. A series of mesoscale 
pressure systems developed along this boundary and moved 
rapidly northeastward. Most of the reported tornadoes were 
associated with these mesolows. The tornado events are in- 
dicated on the surface charts with their reported times of 
occurrence. During the afternoon mesoscale pressure pertur- 
bations of -.10 to -.20 inches (4-8 mb) existed along the 
leading edge of the severe convection. By middle to late 
afternoon the mesohigh pressure area to the rear of the major 
squall-line became stronger with surface pressure perturbations 
in excess of +0.04 inches (>. 2 .mb), 

By 1800 GMT (Fig. 24b) several tornadoes had occurred 
with the mesolow center that was moving northeastward over 
northwestern Alabama. Strong and gusty southerly flow was 
indicated in Mississippi, Alabama, and Georgia to the south 
and east of the squall-line. All of the significant convective 
storms continued to occur well in advance of the synoptic 
scale cold front. At 1800 GMT (Fig. 24c) a major mesolow 
pressure area was present in eastern Mississippi with pressure 
perturbations as great as -0.20 inches indicated. The first 
mesoscale low pressure center had moved to northeastern 
Alabama with two tornadoes reported in this region between 
1800 and 1900 GMT. Pressures were reported to be falling 
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Fig. 24a. Surface mesoanalysis for 1200 GMT, 4 April 1977. 

Pressure analysis (of altimeter settings), fronts 
trough lines, and thunderstorm outflow boundaries 
are solid black lines. Analyses of mesoscale 
pressure perturbations (dashed lines) are shown 
with areas of £ -0.08 inches shaded and areas of 
_> +0.04 cross-hatched. Tornado events logged at 
the National Severe Storms Forecast Center are 
indicated (shaded triangles) along with the GMT 
time of occurrence. 


rapidly at many stations in Alabama, northern Florida and 
western Georgia. There was a marked temperature contrast 
across the squall-line with mid-eighties reported over south 
eastern Alabama compared with upper-sixties temperatures 
over northeastern Alabama. A concurrent GOES visible image 





Surface mesoanalysis for 1500 GMT, 4 April 1977 
Details as in Fig. 24a. 



Surface mesoanalysis for 1800 GMT, 4 April 1977 
Details as in Fig. 24a. 
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Fig. 24d. Full resolution GOES visible image for 1800 GMT, 

4 April 1977 . 

is shown in Fig. 24d. The mesolow in eastern Mississippi was 
located on the southern edge of a large complex of storms, 
however, the low pressure perturbations over northern Alabama 
appeared to be located in a region of suppressed convection 
beneath thin cirrus anvil cloud. 
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Around 1800 GMT two separate mesolow pressure centers 
developed within the low pressure region over eastern Mississippi. 
Figures 24e, f, and g indicate that both of these mesoscale 
pressure systems moved northeastward at 60 to 65 kt ! Indeed, 
Miller (1978) reported that radar cell motions for the severe 
storms on this afternoon averaged about 63 kt . The leading 
center became the second mesolow pressure system to cross 
north-central Alabama. This center was located near Tuscaloosa, 
Alabama, at 1900 GMT; just northeast of Birmingham by 2000 GMT; 
and slightly southwest of Rome, Georgia, by 2100 GMT. 

(Tuscaloosa is located approximately 100 km west-southwest of 
Birmingham and Rome is about 100 km northwest of Atlanta.) 

The trailing system moved along a similar path about one hour 
later. 

Both mesolows were tornadic at 2100 GMT (Fig. 24g) with 
one located near Birmingham and the other just southwest of 
Rome. The region of rapidly falling surface pressures was 
located to the east and northeast of these mesosystems. A 
strong temperature gradient was associated with the squall- 
line and very strong and gusty southerly flow was indicated 
in the warm air mass. A concurrent satellite photograph of 
the thunderstorm activity is shown in Figure 24h. The mesolows 
were associated with intense storms characterized by multiple 
overshooting tops. The high tops seem to indicate a double 
wave structure in Alabama and Georgia with one wave crest 
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Fig. 24e. Surface mesoanalysis for 1900 GMT, 4 April 1977. Details 
(except for pressure perturbations) as in Fig. 24a. 



Fig. 24f . Surface mesoanalysis for 2000 GMT, 4 April 1977. 

Details (except for pressure perturbations) as in 
Fig. 24a. 


58 




located just west of Birmingham and the other just southwest 
of Rome, Georgia. 

By 0000 GMT (Fig. 24i) the major squall-line had moved 
across most of northern and western Alabama and northern Georgia. 
The second and third mesolows that had crossed northern Alabama 
apparently merged as they continued to move northeastward to 
western North Carolina. The magnitude of the low pressure per- 
turbation had decreased markedly. Several final tornadoes were 
reported as this mesosystem moved rapidly northeastward. 


i 
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Fig. 24g. Surface mesoanalysis for 2100 GMT, 4 April 1977. 
Details as in Fig. 24a. 
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Fig. 24h. Full resolution GOES visible image for 2100 GMT 
4 April 1977 . 


Time series showing the surface observations at Birmingham 
and Rome are presented in Figs. 25 and 26. The Birmingham 
figure shows that the gusty southerly winds shifted to the 
southwest and decreased in speed as the first thunderstorm 
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activity moved across the station from 1700 to 1800 GMT. The 
pressure then fell rapidly during the next two hours as the 
leading (second) mesolow moved out of western Alabama and 



Fig. 24i. Surface mesoanalysis for 0000 GMT, 5 April 1977. 
Details as in Fig. 24a. 
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approached the station (See Figs. 24e and f ) . During this time 
the winds backed to a southerly direction and became very strong 
just before the mesolow passed the station. The pressure rose 
rapidly as thunderstorms with the mesosystem affected the station 
and then fell rapidly as the trailing (third) mesolow approached 
from the southwest. The F5 tornado was embedded in this third 
mesosystem. After it passed at about 2100 GMT surface winds 
remained west to northwest and pressures gradually rose as 
continuing thunderstorm activity affected the Birmingham area. 
Notice that during the entire period of extreme pressure fluctu- 
ations the surface temperatures and dewpoints remained relatively 
constant (this is similar to results already shown for the Omaha 
and Neosho tornadoes). 

At Rome, thunderstorms moved over the station around 1900 
GMT in conjunction with the first mesolow and the wind shifted 
to the southwest and the temperatures dropped sharply as the 
dewpoint increased. Surface pressures leveled off for about 
an hour and then fell rapidly as a second mesolow pressure 
system approached. The second center passed the station shortly 
after 2100 GMT and the low pressure remark of 28.74 inches at 
2110 indicates that a tornado cyclone must have passed very 
near the observation site. The surface winds remained quite 
light during this event. The third mesolow passed between 
2200 and 2300 GMT, after spawning the Birmingham tornado. 

The fact that the presence of a mesolow pressure system 
increases the likelihood of severe convective storms has been 
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Inches 



Fig. 25. Time series of surface observations for Birmingham, 
Alabama, on the afternoon of 4 April 1977. Winds 
are in knots with a full barb = 10 kt . 
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known for some time (Magor, 1958 and 1971) . Hoxit etal. 


(1976) suggested that some mesolows may be a hydrostatic 
reflection of middle and upper-level subsidence warming around 
and downwind of thunderstorms . Earlier studies by Williams 
(1957) and Feteris (1961) supported the concept of compensating 
mesoscale circulations in the thunderstorm near environment 
and Fritsch (1975b) and Hoxit et al . (1976) both discussed types 
of feedback processes that might result. 

A number of analyses of upper-air features during the 
4 April 1977, event are presented in Figs. 27, 28, and 29. 
Isotherm analyses are for layer mean temperatures. Since the 
thunderstorm activity occurred well ahead of the synoptic scale 
front and its associated upper-level baroclinic zone, the changes 
in temperature that were found in the region of thunderstorms 
may be partially attributed to convective processes. 

Figures 27a-c present analyses of lower tropospheric 
conditions for the period 1200 to 0000 GMT. Winds and heights 
are plotted for the 850 mb level and 850 mb streamlines are 
indicated. Mean temperatures are plotted, and isotherms drawn, 
for the surface to 700 mb layer. During the period, the strong 
baroclinic zone associated with the synoptic cold front advances 
slowly eastward across the Mississippi Valley. However, a cool 
pocket and a strong baroclinic zone develop in association with 
the squall-line. There is no way that advective processes 
could totally account for the cooling observed over central 
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Fig. 27a. Lower tropospheric conditions for 1200 GMT, 4 April 
1977. Heights, winds, and streamlines are for 850 
mb. Isotherms are shown for the mean surface to 
700 mb temperature (plotted below the height). 


Alabama between 1800 and 0000 GMT (compare Figs. 27b and c) . 

The eastern baroclinic zone is a mesoscale feature that has 
been generated by the storm activity. Low-level cooling is 
the primary process responsible for the development of thunder- 
storm mesohigh pressure systems. However, the hydrostatic 
surface pressure is a reflection of the temperature profile 
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Fig. 27b. Lower tropospheric conditions for 1800 GMT, 

4 April 1977. Details as in Fig. 27a. 

through the depth of the atmosphere and a mesohigh can form 
in the cool, lower tropospheric air only if the cooling is not 
offset by upper-level warming. 

Middle-level conditions and changes are illustrated in 
Figs. 28a-c. These charts show heights, winds, and streamlines 
for the 500 mb level and thermal analyses for mean temperature 
in the 700 to 250 mb layer. During the analysis period strong, 
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Fig. 27c. Lower tropospheric conditions for 0000 GMT, 
5 April 1977. Details as in Fig. 27a. 
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middle-level warming occurs in the region of the severe squall- 
line. Mean temperatures in the layer increase from 3-5°C 
during the six hours from 1200 to 1800 GMT over portions of 
northern Mississippi, northern Alabama, and central Tennessee. 
Part of the increase is due to warm advection; however, the 
tremendous amplification of the thermal ridge during the twelve 
hour period suggests that the warming is partly due to the 
large area of convective storms. This warming could be caused 
by two primary processes. The first is direct warming due to 
persistent moist adiabatic ascent over a large area. Secondly, 
the warming could also be produced by compensating subsidence 
within the environment surrounding the thunderstorm complexes. 

Upper tropospheric and lower stratospheric conditions 
and changes are illustrated in Figs. 29a~c. These charts show 
heights, winds, and streamlines for the 200 mb level and thermal 
analyses for mean temperature in the 225 to 175 mb layer. 
Although warm advection was indicated over Alabama and eastern 
Tennessee (see Figs. 29a and b) a large area of significant 
upper tropospheric cooling occurred in the region of convective 
activity (see Fig. 29c). Mesoscale lifting and cooling at the 
cloud-top level was probably the mechanism responsible for the 
observed changes. 

Thermal changes similar to those illustrated in Figs. 27, 

28 and 29 have been obtained by Fritsch (1978) in three- 
dimensional numerical simulations of squall-line convection. 
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Fig. 28a. Middle tropospheric conditions for 1200 GMT, 

4 April 1977. Heights, winds, and streamlines 
are for 500 mb. Isotherms are shown for the 
mean 700 to 250 mb temperature (plotted below 
the height) . 
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Fig. 28b. Middle tropospheric conditions for 1800- GMT, 
4 April 1977. Details as in Fig. 28a. 
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Fig. 28c. Middle tropospheric conditions for 0000 GMT, 
5 April 1977. Details as in Fig. 28a. 
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These figures illustrate that a large region of organized 
convective storms can affect and modify their near environment 
on scales large enough to be detected in the synoptic sounding 
network. Convective modification can be easily detected at 
the surface where observations are fairly dense and the convec- 
tive changes dramatic. At middle and high levels, however, 
convective mesoscale features are more difficult to detect 
because of large data separations and the superposition of 
mesoscale features upon strong background flow fields. 

The large area of higher pressure from central Alabama 
into eastern Tennessee (see Figs. 24g and 24i) appears to be 
directly related to the area of lower tropospheric cooling 
evident on Fig. 27c. The exact mechanisms responsible for 
the series of mesolows, however, are not resolved by the 
coarse resolution, upper-air data. The April 4th mesolows 
appeared to be closely linked to active storm complexes; 
whereas, the April 24th mesolow (considered earlier) developed 
in a warm lower tropospheric region prior to storm development. 

The April 4th mesolows were probably a result of a combination 
of mesoscale effects that would include: low-level warm advection, 
strong moist adiabatic warming, and upper-level subsidence 
and warming . 

b. Upper- tropospheric flow fields oh 24 and 25 April 1975 
Maddox (1979a) has examined the evolution of jet stream 
level winds using the NASA AVE IV data. His findings provide 
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Fig. 29a. Upper tropospheric/lower stratospheric conditions 
for 1200 GMT, 4 April 1977. Heights, winds, and 
streamlines are for 200 mb. Isotherms are shown 
for the mean 225 to 175 mb temperature (plotted 
below the height) • 
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Fig. 29b. Upper tropospheric/ lower stratospheric conditions 

for 1800 GMT, 4 April 1977. Details as in Fig. 29a. 
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Fig. 29c. Upper tropospheric/lower stratospheric conditions 

for 0000 GMT, 5 April 1977. Details as in Fig. 29a. 
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excellent examples of convective feedbacks that affect the larger 
scale environment and they are briefly summarized in this section. 

Figure 30a shows an analysis of jetstream winds (defined 
as the maximum reported wind above the 350 mb level) for 0000 
GMT, 24 April 1975. Isotachs and jet axes are shown along With 
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Fig. 30a. Jetstream winds (maximum wind above 350 mb) for 
0000 GMT, 24 April 1975. Isotachs (kt) and jet 
axes are shown along with radar echoes from the 
2335 GMT radar summary chart (from Maddox, 
1979a). 
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radar echoes from the radar summary facsimile chart. The 
analysis shows a diffuse polar jetstream from western Kansas 
to New England with maximum speeds on the order of 100 kt . 

The subtropical jet, also with maximum speeds of about 100 kt , 
curves from southwest Texas across Georgia. The 0000 GMT 
digitized IR satellite image -(Fig. 30b) shows that intense 
storms are occurring over northern Missouri and southern Iowa 
(radar tops 50-60,000 ft.). The jetstream chart for 0600 GMT 
24 April 1975 (Fig. 31a) shows that tremendous changes have 
occurred during a period of only 6 hours. The subtropical jet 
is poorly defined as it crosses the Gulf Coast states. However, 
a significant jet, with maximum speeds exceeding 130 kt , curves 
from eastern Nebraska to the East Coast. Several secondary jets 
are also indicated. The strong flow has developed along the 
north edge of a region of strong thunderstorms (radar tops 
40-50,000 ft. from Michigan to Kansas). The flow splits and 
diverges around the western end of the convective system as 
if a mesohigh pressure perturbation exists over the storm 
complex. The IR satellite image for 0600 GMT (Fig. 31b) shows 
that the strong jetstream is positioned along the northern 
edge of the cirrus shield associated with the thunderstorm 
complex . 

Figure 32a shows the jetstream chart 18 hours later at 
0000 GMT on the 25th of April, The subtropical jet was once 
again well defined with speeds in excess of 125 kt indicated 
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Fig. 30b. Digitized IR satellite image for 0000 GMT, 
24 April 1975. 
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Fig. 31a. Jetstream winds for 0600 GMT and radar echoes 
for 0535 GMT, 24 April 1975. Details as in 
Fig. 30a (from Maddox, 1979a). 
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Fig. 31b. Digitized IR satellite image for 0600 GMT, 

24 April 1975. 

over central Texas. A diffuse polar jet stretched over the 
northern U.S. with speeds of about 100 kt. The shower activity 
and slightly higher wind speeds over the East Coast were the 
remains of the convective system and jet shown in Fig. 31a. 
Intense thunderstorms were developing over western Missouri 
and eastern Oklahoma at the time of the figure. The satellite 
IR photograph shown in Fig. 32b shows these developing storms. 
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Fig. 32a. Jetstream winds for 0000 GMT and radar echoes 

for 2335 GMT, 25 and 24 April 1975. Details as 
in Fig. 30a (from Maddox, 1979a). 


The wind field at the jetstream level six hours later 
(Fig. 33a) shows that tremendous changes have again occurred 
in only six hours. By 0600 GMT a large, intense thunderstorm 
complex (radar tops 40-55,000 ft.) extended from Kentucky 
westward to eastern Oklahoma. A jetstream with maximum speeds 


in excess of 135 kt curled anticyclonically around the northern 
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Fig. 32b. IR satellite image for 0000 GMT, 25 April 1975. 


edge of the thunderstorm system and the extreme diffluence over 
the system again suggests the presence of a mesohigh pressure 
perturbation aloft. The concurrent IR satellite image (Fig. 
33b) shows the intense storm system and indicates that the 
jetstream maximum is again located along the northern- edge of 
the extensive cirrus shield. 
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Fig. 33a. Jetstream winds for 0600 GMT and radar echoes 
for 0535 GMT, 25 April 1975. Details as in 
Fig. 30a (from Maddox, 1979a). 

Maddox (1979b) has documented the existence of a large 
mesohigh circulation at 200 mb at the top of a similar con- 
vective system using an objective analysis scheme to isolate 
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mesoscale features. Fritsch et al.(1979) have shown results 


from numerical simulations of squall-line type convection that 
are similar to the features observed on 4 April 1977 and 24/25 
April 1975. Ninomiya (1971a and 1971b) used satellite photo- 
graphs to document similar upper-level wind perturbations. 

All of these examples show that when convective storms organize 
into large mesosystems the convectively driven circulations 
affect and modify the large scale flow fields on scales that 
may, on occasion, be detected in the synoptic upper-air data. 



Fig. 33b. Digitized IR satellite image for 0600 GMT, 
25 April 1975. 
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5. 


SUMMARY 


Results from several detailed case studies of severe storm 
events have been presented to document various ways in which 
convective storms interact with their environment. A physical 
model of boundary layer wind fields (both in the horizontal 
and vertical) was developed to explain why thunderstorms 
intensify as they move along or across thermal boundaries. 
Detailed time series data, along with three specific case 
studies, substantiated many aspects of the proposed model. 

Winds veer slowly with height on the cool side of such boundaries 
while winds veer rapidly through the subcloud layer in the hot, 
moist airmass. If mean subcloud winds are considered, meso-g 
scale convergence and cyclonic vorticity tend to obtain maximum 
values in a narrow "mixing zone" along the thermal boundary. 

The 1977 examples demonstrate that this local intensification 
effect is important when large scale conditions do not favor 
family outbreaks of severe storms. In these situations "mini- 
outbreaks" of strong tornadoes occur along thermal boundaries. 

Objective analyses of surface kinematic fields were mapped 
onto GOES satellite imagery. It was demonstrated that meso-g 
scale areas where conditions were favorable for significant 
tornadic storms could be identified and that the satellite 
imagery could be used to monitor storm development in, or 
movement into, these dynamically favorable regions. Severe 
storm nowcast techniques based on these results should be tested 
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and studied for a number of severe storm situations. An 
interactive data processing and analysis facility would be 
desirable for use in further studies and absolutely necessary 
for any real-time tests of applications. 

The final examples presented demonstrated that once severe 
thunderstorms organize into meso-a scale storm complexes the 
associated mesoscale circulations act to modify larger scale 
features. Mesoscale pressure systems, both at the surface and 
aloft, may develop. Once formed, these systems may act to 
enhance and maintain the convective storms. Many of these 
important features and circulations exist aloft since most of 
the examples showed large mesoscale pressure fluctuations 
occurring without corresponding fluctuations in surface tempera- 
ture and dewpoint. The cases that show pronounced diffluence 
and jetstream maxima developing in response to large storm 
complexes are important. These patterns are usually thought 
of as precursor features that favor storm development. Con- 
vective feedbacks, such as those illustrated, cannot occur in 
current operational numerical models because of the simple 
convective parameterizations that are used. 
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APPENDIX I 


The Barnes Analysis Technique 

The low-pass filter, exponential weight function objective 
analysis technique developed by Barnes (1964) was used to 
analyze surface data for the VTP computations presented in 
Sec. 3. The basic objective analysis technique follows Barnes’ 
(1973) revision of his earlier scheme. Two dimensional meteor- 
ological data are interpolated to a cartesian grid and simul- 
taneously smoothed. The first estimate of the gridded values 
of some parameter f is given by: 

N N 

f (i,j) = w f / w (Al) 

o n n *—> n 

n=l n=l 

where the weighting function, w n , is: 

w = exp (-r^/4 k) . (A2) 

n n 

The weight function parameter k is a constant arbitrarily chosen 
to fit the specific application being considered. The distance 
from the (i,j) grid point to the location of the data values 
f is denoted by r and N is the number of data observations 
allowed to influence the grid value. 

Barnes' (1973) modification requires a correction pass 
through the first estimate field to achieve a desired response 
at specific wavelengths. The value of the weight function 
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parameter is reduced on the correction pass resulting in 
decreased computation time compared to the earlier use of 
successive iterations. The final grid point values are 
given by: 

N N 

= f Q (i 5 j) +£ w n D n / £ W A 5 (A3) 

n=l n= 1 

where the modified weight function, w’ , is: 

2 

w^ = exp (-r /4g k) , 0<g<l; and (A4) 

D n = f n ( x , y ) - f Q ( x , y ) . (A5) 

Thus D n is the difference between an observed data value 
and the first estimate value at the same point. The first 
estimate value must be interpolated using f o (i,j) at the four 
nearest grid points. The final data field has been interpolated 
and filtered. Barnes has shown that the filter response is 
given, as a function of wavelength A, by: 

R = R (1 + R g_1 - R g ), where: (A6) 

o o o 

R o = exp (- 7T 2 4k/A 2 ). ( A7 ) 

Advantages of this exponential weight, objective analysis 
technique listed by Barnes (1973) include: 

lv The weight factors k and g are chosen prior to the 
analysis so that pattern scales supportable by the 
data distribution will be revealed to a known response 
amplitude . 
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2. Since (A2) approaches zero asymptotically, the 
influence of data can be extended any distance without 
changing the weight function, and therefore, the 
response characteristics. 

3. Small scale noise is adequately suppressed so that 
further smoothing with numerical filters is not 
necessary . 

4. The desired pattern resolution can be achieved in 
only one iteration. 

Barnes has also listed a number of practical guidelines 
to be used in specifying the values of k and g to fit particular 
applications. The density and distribution of the data field 
being analyzed is the most critical factor affecting the analysis. 
The response function (computed from equation A6) for the 
analysis runs used in this study is shown in Fig. A1 . 
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